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In the presence of a quark chemical potential, a magnetic field induces an axial current 

^—>. ■ in the direction of the magnetic field. We compute this current in the Sakai-Sugimoto model, 

fvj ' a holographic model which, in a certain limit, is dual to large- A'^c QCD. We also compute 

the analogous vector current, for which an axial chemical potential is formally introduced. 

r^ , This vector current can potentially be observed via charge separation in heavy-ion collisions. 

After implementing the correct axial anomaly in the Sakai-Sugimoto model we find an axial 

current in accordance with previous studies and a vanishing vector current, in apparent 

contrast to previous weak-coupling calculations. 

§1. Introduction 
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^ , The Sakai-Sugimoto model^^'^^ is a realization of the gauge/gravity duahty^^ 

QJ ' which provides a top-down approach to a holographic model of QCD. While at 

f^, ' present no gravity dual to QCD is known, the Sakai-Sugimoto model exhibits some 

of the most important properties of QCD, most notably confinement and chiral 

symmetry breaking. Confinement is realized within a background geometry of Nc D4- 

branes, extending in the 3-1-1 dimensions of the field theory and an extra dimension 

Q^ i x^.^^ This extra dimension is compactified on a circle with radius M^^^, where the 

'^ ' Kaluza-Klein mass Mkk is a parameter of the model and sets the scale for the mass of 

. ■ unwanted adjoint scalars and fermions. Two solutions for the background geometry 

f~^ ' between them at a critical temperature Tc = Mkk/CSvt). 

Fundamental chiral (and massless) fermions are described by introducing Nf D8- 
and D8-branes (here and in most applications treated as probe branes), separated 
in the extra dimension. ^^ The gauge symmetries on these branes account for the 
global chiral symmetry SU{Nf)L x SU{Nf)ji in the associated field theory. Two 
5^ \ qualitatively different embeddings of the flavor branes are interpreted as chirally 

symmetric (D8- and D8-branes separated) and chirally broken (D8- and D8-branes 
connected) phases. In the version of the model used here, chiral symmetry breaking 
and confinement are equivalent, i.e., the ground state is chirally broken if and only if 
it is confined (by choosing a sufficiently small distance between D8- and D8-branes 
in the extra dimension one finds a chirally broken, deconfined phase; however, from 
the geometry of the model it is obvious that there cannot be a chirally restored, 
confined phase). 

Here we are interested in the evaluation of the model in the presence of a quark 
chemical potential and a magnetic field. It has been shown that a magnetic field 
gives rise to baryon number even in the chirally broken phase^' ^-^ (the resulting 
baryon density being homogeneous, in contrast to instanton-like baryons in the Sakai- 

typeset using VTpT^.ds (Ver.0.9) 
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Sugimoto niodel^^). In our context, two physical situations are of special interest. 
First, large magnetic fields can be found in the interior of compact stars (surface 
magnetic fields of magnetars being of the order of i? < 10^^ G), where also the quark 
chemical potential is large, n ~ (400 — 500) MeV. In this case, we may expect an 
axial current in the direction of the magnetic field. ^^' ^^' Second, large magnetic fields 
can be temporarily created in noncentral heavy-ion collisions, reaching values of the 
order of B r^ 10^^ G. In the latter case, the interplay of the magnetic field with 
a nonvanishing chirality N^, i.e., a nonzero difference between the number of left- 
and right-handed fermions, may lead to the so-called chiral magnetic effect. ^^''^^' 
In the chiral magnetic effect, an electric current is generated in the direction of 
the magnetic field, which possibly is responsible for the charge separation that has 
been observed at the Relativistic Heavy-Ion Collider (RHIC).^^' The nonzero N^ is 
provided by the QCD axial anomaly which relates certain gluon configurations to 
a change in N^,, and thus to a nonzero current on an event-by-event basis. Here 
and in the original weak-coupling calculations, the gluon field fluctuations are not 
taken into account explicitly, a nonzero Nr, is rather described by introducing an 
axial chemical potential n^. This description has to be taken with some care since 
Nr, is not a conserved quantity and thus strictly speaking there is no such thing as 
an axial chemical potential in the thermodynamic sense. 

In summary, we shall be interested in the axial current in the presence of a 
vector chemical potential and the vector current in the presence of an axial chemical 
potential. More details about the following calculations can be found in Ref.^^' 

§2. Currents and anomalies 

We start from the (Euclidean) action of the Sakai-Sugimoto model for one quark 
flavor,^) 

S = Sym + Scs , (2-1) 



with Yang-Mills (YM) and Chern-Simons (CS) contributions 



d'^x / dz 



N, 



— oo 
oo 



^ 2M? ^ 



(2-2a) 



5cs = ^ / d^x I dz A^F^.Fp^e^'^P- , (2-2b) 



— CO 



with Greek indices running over fi,iy, ... = 0,1, 2, 3. Our convention for the epsilon 
tensor is eoi23 = +1- We have used the metric functions 

kiz) = l + z\ hiz) = (1 + z^)-^/\ (2-3) 

and the dimensionless constant 

(2.4) 



2167r3 ' 

where A is the 't Hooft coupling. The integration is over space-time (r, x) and the 
holographic (dimensionless) coordinate z. In the confined phase, z extends from the 
left-handed boundary {z = +oo) over the tip of the cigar-shaped {x^,z) subspace 
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(z = 0) to the right-handed boundary (z = — oo). (In the deconfined phase, the 
integral over connected D8 and D8 branes is replaced by two separate integrals over 
disconnected D8 and D8 branes; in this section we restrict ourselves to the confined 
phase, all arguments are analogous in the deconfined phase.) 
The equations of motion are 

kM^k d,[kiz)F'^] + '^h{z)d,F''^' = -^F^.Fp^e^'-P- , (2-5a) 

kM|k d^[k{z)F^'^] = ^F^,F,,e>^^P- , (2-5b) 

where the second equation is obtained from varying A^ (which, in our gauge choice, 
has already been set to zero in the action (2-2)). 

Next we introduce the chiral currents. They are defined through the variation 
of the on-shell action with respect to the boundary values, 

(2-6) 
where the first (second) term is the YM (CS) contribution. It is only the YM part 
of the current, i.e., the first term in eq. (2-6), which appears in the asymptotic 
expansion of the gauge fields:^^ from the definition (2-6) and with k{z) = 1 + z^ we 
find 

j^/^ 1 / 1 \ 

^„(.,3)=A.(-.,. = ±co)±^- + o(^). (2.7) 

One can also confirm this relation from our explicit results below. 

The divergence of the currents (2-6) can be easily computed with the help of the 
equation of motion for A^ (2-5b). With the left- and right-handed field strengths 
F^xl (x) = F^,y{x,z = ±oo), the corresponding dual field strength tensors i^^T^ = 

i FpJ^e'^'P'', the vector and axial currents J'' = Jl^ + J^, J^ ^ Jr - Jl, and the 
vector and axial field strengths introduced as F^ = F^ + Fp^, F^j^ = F^ — F^, we 
obtain the vector and axial anomalies 

d,Jt = ^ (<.^r + Fi^Fif) ■ (2-8b) 

The coefficients on the right-hand side (which receive contributions from both the 
YM and CS parts of the currents) are in accordance with the standard field-theoretic 
results for Nc chiral fermionic degrees of freedom coupled to left and right chiral 
gauge fields. ^^^ The above form of the anomaly, which is symmetric in vector and 
axial- vector gauge fields, is called consistent anomaly. If left- and right-handed Weyl 
spinors are treated separately, this form of the anomaly arises unambiguously. In 
QED, however, we must require that the vector current be strictly conserved, even 
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in the presence of axial field strengths. As was first discussed by Bardeen/^' this can 
be achieved by the introduction of a counterterm that mixes left- and right-handed 
gauge fields. Having even parity, Bardeen's counterterm is uniquely given by^^^ 

AS = cjd'xiAf^A^F;^^ + ^^^^F^^Je'^'^^- , (2-9) 

where c is a constant determined by requiring a strictly conserved vector current. 
Because this expression can be naturally written as a (metric-independent) integral 
over a hypersurface at |z| = yl — )■ oo with left- and right-handed fields concen- 
trated at the respective brane locations, AS can actually be interpreted as a (finite) 
counterterm in holographic renormalization. In particular, it does not change the 
equations of motion. 

The contribution of Bardeen's counterterm to the chiral currents is 

AJj",^ = ^c (a^/^FJ^J^ - A^'^FJ^J^ + 2A^'^F^^^) e^"^'^ . (2-10) 



L/R T^y^'^u ^ pa ^'-v ^ pa " ^^'-v ^ pa 

Denoting the renormalized left- and right-handed currents by J^in = •^L/p + ^-^r/i?' 
and J^ = J^^jf^, J^ = J^- J^, we find that c = 7Vc/(487r2) is required to obtain 
the covariant anomaly 



df,J'' = 0, (2- 11a 






The prefactor in front of the first term in the axial anomaly now has changed from 
A^c/(247r2) in eq. (2-8b) to Nc/{8tt'^), which is the well-known result for the Adler- 
Bell-Jackiw anomaly for QED^^-''^^-' and which is essential for getting the correct 
pion decay rate vr" — )■ 27. The necessity of adding the counterterm (2-9) to the 
Sakai-Sugimoto model is in fact completely analogous to the well-known procedure 
in chiral models where a Wess-Zumino-Witten term accounts for the anomaly. "^^•' 

In the literature sometimes the coefficient of the subleading term in the asymp- 
totic behavior oi A^(z) and thus the YM part of the current (see eq. (2-7)) is identified 
with the full current. 2*^) '^i) Using this identification, it has also been assumed that 
the equation of motion for A^ (2-5b) represents the anomaly equation. ^2) Indeed, 
from eq. (2-5b) one obtains 

^.^YM,5 = ^ (^.^.^r + ^A^r) , (2- 12b) 

and this does contain the same coefficient in front of F^^Fy'^ as the full covariant 
anomaly (2-11). However, it differs from the latter in the presence of axial gauge 
fields. In particular, the vector current is then not strictly conserved. Even when 
this issue may be ignored, because all axial vector field strengths are set to zero, it 
appears to be questionable to keep only part of the full current (2-6). 
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§3. Solution to the equations of motion 

We introduce a homogeneous magnetic field B in the spatial 3-direction via the 
gauge field component Ai{x2) = —X2B, i.e., Ai is a constant with respect to the 
holographic coordinate z, which is consistent with the equations of motion. Then, 
the only nontrivial gauge field components we need are ^o(^) and ^3(2), where Aq is 
necessary to account for finite chemical potentials, and ^3 becomes nonzero through 
the equations of motion. (With only a little more eff'ort, the equations of motion can 
also be solved in the presence of an additional electric field. ^^-^ This is for instance 
useful to check the axial anomaly explicitly, but for simplicity we restrict ourselves 
to a magnetic field which is the physically relevant field in our context.) 

3.1. Chirally broken phase 

Within this ansatz, the equations of motion read 

d,{kd,Ao) = 2f3d,A3 , d,{kd,As) = 2(3d,Ao , (3-1) 

with the dimensionless magnetic field /3 = aB/M^^ where a = 27tt/{2X). The 
general solutions to Eqs. (3-1) are 



sinh(2/3 arctan z) 
Ao{z) = fi- fi^- 



^3{z) = -/U5 



sinh (3-ir 
cosh(2/3 arctan z) 



cosh(2/3 arctan z) 

coth pvr 



sinh /Stt 



sinh /3-iT 
coth fSir 



(3-2a) 



sinh(2/3 arctan z) /in 

J r-— 3 . (3-2b) 

smh pn 



Here we have identified the boundary values of the temporal components of the gauge 
fields with the chemical potentials, Ao{z = ±00) = fiL,R, and the vector and axial 
chemical potentials are fi = {^r + fJ-L)/"^ and ^5 = {fiR — ^l)/2. The boundary 
values of the spatial components of the gauge fields are A^lz = ±00) = =FJ- The 
axial supercurrent j, which has to be determined from minimizing the free energy, 
describes a rotation of the chiral condensate in the space spanned by the scalar and 
pseudoscalar mesons^) and thus may be related to the recently discussed quarkyonic 
chiral spirals. ^^^ 

3.2. Chirally restored phase 

In the case of disconnected D8 and D8 branes, i.e., in the deconfined, chirally 
restored phase, we have two sets of equations of motion, one for the left-handed fields 
A^{z) and one for the right-handed fields A^{z) with z G [0, 00], 



'"/^^ - -f a.Aj/^ , 9.(^35.4/^) = ±f ^ ^^/^ 



d^ikodX^ ) = ±^^.4^ . d^ik^d^A^^'') = i^d^^"" . (3-3) 



Here 6 = 27rT/MKK is the dimensionless temperature, and ko{z) = - — ^ — , ks{z) = 
z{l + z'^Y''^ are diff'erent metric functions for temporal and spatial components of 
the gauge fields. The general solution to Eqs. (3-3) is 



4/^(z) = (^TM5) 



p{z) q{z) 

qo 



(3-4a) 
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^/^(z) 
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kodzP {1 + kodz 

Qo 



where po = p{0), qo = q{0), and 

Vl - 16/3'2 + 1 Vl - 16/3'2 - 1 1 1 



p(z) 
q{z) 



2-f'l 



2'l + z2 



1 



\/rT^ 



2F1 



^\ - 16/3'2 - 1 ^1 - 16/3'2 + 1 3 1 



2'l + z2 



(3-4b) 

(3-5a) 
(3-5b) 



with the abbreviation /?' = P/O''^ . Again, the boundary values of the temporal com- 
ponents of the gauge fields are identified with the chemical potentials, ^q (00) = 
fJ^L/n- Following Ref.,2'^) these components vanish at the horizon, A^^' (0) = (see 
however Refs.^^^'^^'). The spatial components vanish at the holographic boundary, 
but acquire a finite value at the horizon, A^ (0) 






§4. Results and discussion 

To compute the axial and vector currents in the direction of the magnetic field 
we use the solutions to the equations of motion, Eqs. (3-2) and (3-4), and insert them 
into the chiral currents, which are defined through Eq. (2-6) and the contribution 
from Bardeen's counterterm (2-10). The results are summarized in Table I. 

Let us first discuss the full renormalized currents J = i7ym + iTcs + ^^ which 
have been obtained under the constraint of the correct (covariant) QED anomaly. 
For the axial current we find that the counterterm A J' exactly cancels the CS part. 
In the chirally symmetric phase, this yields exactly the expected topological result 
iTs = ^^Y^- ^^ particular, the dimensionful parameter of the model, MkK) 

drops out of this universal result. In the chirally broken phase, the axial current is 
smaller and given by a complicated function of B. 

The most striking of our results is that for both phases the renormalized vector 
current is zero for all magnetic fields. One might have expected that, in the decon- 
fined phase, this current should reproduce the known result for the chiral magnetic 
effect, ^2) J' = t^5 ^ e ^ because this result can be derived from the anomaly only^'^) 
and we have explicitly made sure that our model reproduces the correct anomaly. 
The question is thus whether our result can indeed be interpreted as the current 
responsible for the chiral magnetic effect. We shall comment on possible issues re- 
lated to this question below. Our result of a vanishing current in the confined phase, 
however, seems less puzzling. The usual explanation of the chiral magnetic effect, 
using a quasiparticle picture (which is not guaranteed to hold in our strong-coupling 
approach), relies on individual, electrically charged, massless quarks which move in 
different directions according to their chirality. A suppression of the effect may thus 
indeed be expected in the confined, chirally broken phase. ^^^'^^-' 

In Table I we also show the results which are obtained from the YM part only. In 
the case of the axial currents, there is no difference to the renormalized currents. For 
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Jym 


Jym + Jcs 


Jym + Jcs + A J 




anomaly 


"semi-covariant" : 


consistent: 


covariant: 


d.Jt/^ 


3FvFv+3FaFa 


FvFv + FaFa 


ZFvFv + FaFa 


9.J'/^ 


GFvFa 


2FvFa 





T>Tc 


j^r-§^ 


1 


2 
3 


1 


Jl^^i^ 


1 


2 
3 





T <Tc 


J,l^ 


/3 coth/37r 
2p(/3) 


/9coth/37r 1 
2p(,9) 3 


;3coth,3ir 
2p{;3) 


JI^^BN^ 


1 


2 
3 






Table I. Different (parts of the) axial and vector currents Jr> and J - normalized to the weak- 
coupling results®'' ^^^ - in the direction of the magnetic field B in the deconfined, chirally sym- 
metric phase (r > Tc) and in the confined, chirally broken phase (T < Tc). We have abbreviated 



p(/3) = /3coth/37r -I- ^'^.^ „ with the dimensionless magnetic field /? — aB/My. 



the vector currents, however, we now recover the weak-couphng (deconfined) result 
for both confined and deconfined phases. Therefore, now the deconfined current is 
as expected while the lack of suppression in the confined phase comes as a surprise. 

There are several problems in the current approach. First, upon computing the 
free energy explicitly and then taking the derivative with respect to the appropriate 
source, the currents turn out to be different from the straightforward definition via 
the gauge/gravity correspondence (used for the results in Table I). This disturbing 
discrepancy arises for nonvanishing background magnetic fields and can be attributed 
to boundary terms at spatial infinity. A previously suggested fix of this problem by 
modifying the action^' seems to be not acceptable for our purpose because it entirely 
eliminates the axial anomaly from the correspondingly modified currents. ^^' Second, 
the introduction of a /is may be problematic. Since a chemical potential must be 
associated to a conserved charge, it has been argued that a physically meaningful 
vector current can be obtained only after introducing a conserved (not anomalous) 
charge A^5.^^)''^^) This, on the other hand, corresponds to a gauge variant axial 
charge density, which precludes a generalization to inhomogeneous situations. ^^•' 

In a simple holographic model not unlike the Sakai-Sugimoto model, and in a 
setup that corresponds to our deconfined phase, the vector current has recently been 
found to coincide with the weak-coupling result. ^^•' In this study, which uses linear 
response theory (i.e., only deals with an infinitesimally small magnetic field), it has 
been pointed out that it is important to distinguish the source for the quark den- 
sity (= boundary value of the gauge field) from the chemical potential (= potential 
difference between boundary and horizon). With this distinction the sources can be 
set to zero, and consequently the CS contributions to the currents vanish trivially. 
It remains to be seen whether this concept can be applied to the present model, in 
particular to the case of a finite background magnetic field. ^^•' Moreover, the situa- 
tion in the confined phase seems profoundly different, since in the case of connected 
branes it is unclear how to define the chemical potential other than through the 
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boundary value. Maybe this difference is not surprising because it is the deconfined 
phase in which the results of the renormalized currents seem puzzling and the YM 
currents seem to agree with physical expectations. In the confined phase it is the 
other way around. 
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